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1 .0  SUMMARY 

In  an  effort  to  better  understand  the  mechanics  of  wave  propagation 
through  the  near  field  of  an  explosive  seismic  source,  the  real  time  functions 
of  stress  and  strain  have  been  measured  for  a  variety  of  rock  types  using  both 
tensile  and  compressive  loading.  Experiments  were  performed  using  cylindrical 
specimens  of  Boise  sandstone,  Berea  sandstone,  Indiana  limestone,  and  Westerly 
granite.  Peak  stresses  varied  from  4.5  *  10*  Pa  to  2.5  *  10®  Pa,  and  maximum 
longitudinal  strains,  as  measured  on  the  surface  of  the  cylinders,  varied  from 
approximately  1  x  10”®  to  2  *  10“*.  Differences  in  the  details  of  the  experi¬ 
mental  results  exist  among  the  various  rock  types  studied,  but  the  following 
observations  are  generally  applicable: 

1.  At  very  low  amplitudes  the  rocks  behave  almost  linearly;  the 
stress-strain  curves  have  the  same  slope  on  either  side  of  zero 
stress,  and  a  hysteresis  loop  cannot  be  detected  above  the  noise 
level . 

2.  The  first  evidence  of  nonlinearity  appears  at  peak  stresses  of 
approximately  1  *  10®  Pa.  At  this  amplitude  the  slope  of  the 
stress-strain  relationship  is  not  the  same  on  either  side  of 
zero  stress.  In  all  but  two  anomalous  cases,  Indiana  limestone 
and  Westerly  granite  at  relatively  low  peak  stress  levels,  the 
effective  modulus  is  higher  in  compression  than  in  tension. 
Furthermore,  a  small  hysteresis  loop  is  observed,  especially  on 
the  tensile  side  of  zero  stress. 

3.  Hysteresis  loops  continue  to  enlarge  with  increasing  maximum 
stress.  Compresslonal  hysteresis  loops  are  smaller  than 
tensional  loops,  and  nonlinearity  is  obvious. 
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These  results  confirm  that  nonlinear  effects  are  present  in  rocks  at 
very  low  amplitudes  (<105  Pa  stress  and  <1£T5  strain)  both  in  compression  and 
in  tension.  Material  response  in  compression  is  more  elastic  than  in  tension, 
in  that  the  loading-release  curve  follows  a  curved  path,  but  with  relatively 
little  hysteresis.  A  significant  amount  of  hysteresis  is  observed  in  ten¬ 
sion.  This  is  apparently  related  to  intergranular  friction,  which  is  aided  by 
the  opening  of  Intergranular  cracks.  Measurements  using  the  resonating  bar 
technique,  as  reported  previously,  are  excellent  for  detecting  the  first  signs 
of  nonlinearity  in  either  tension  or  extension,  but  cannot  resolve  differences 
in  extensional  behavior  vs  compressional  behavior.  Finally,  it  is  unlikely 
that  rocks  respond  perfectly  linearly  to  high  amplitude  compressive  pulses 
propagating  through  the  near  field  of  an  explosion  at  peak  stresses  greater 
than  approximately  105  Pa.  It  is  more  likely  that  material  response  is 
nonlinear  in  compression,  with  a  small  amount  of  hysteresis  due  to  inter¬ 
granular  friction.  Further  studies  of  rocks  are  planned  while  the  specimens 
are  subjected  to  elevated  effective  stresses  in  order  to  test  this  hypothesis. 
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2.0  INTRODUCTION 

The  mechanics  of  wave  propagation  through  the  near  field  of  an 
explosive  seismic  source  are  not  completely  understood,  particularly  in  the 
peak  strain  amplitude  range  between  10"^  and  10"®.  The  most  complete  data 
sets  which  describe  free-field  wave  propagation  as  a  function  of  scaled 
distance  (distance/yield*/'*)  in  this  strain  amplitude  range  have  been  obtained 
for  explosions  in  polycrystalline  salt.  The  data  base  includes  measurements 
from  (a)  the  5.3  kT  nuclear  event,  SALMON  (Perret,  1967;  Trulio,  1978);  (b) 
the  medium  to  large  scale  chemical  explosions  of  the  COWBOY  experiments,  which 
took  place  in  a  natural  salt  dome  (Trulio,  1978,  1981);  and  (c)  a  series  of 
small  scale  chemical  explosions  in  pressed  polycrystalline  salt  in  the  labor¬ 
atory  (Larson,  1982).  When  all  available  data  are  combined,  wave  propagation 
appears  to  satisfy  cube  root  of  yield  scaling.  That  is,  the  decay  of  peak 
particle  velocity  and  displacement  can  be  defined  as  a  function  of  distance/ 
yield1/3  (Trulio,  1978,  1981;  Larson,  1982).  This  appears  to  be  the  case, 
even  though  available  data  cover  approximately  10  decades  in  yield,  4  decades 
in  peak  particle  velocity,  and  4  decades  in  frequency.  However,  the  data  fit 
a  line  that  decays  much  faster  than  r"1,  which  would  be  the  case  if  the  mater¬ 
ial  behaves  in  a  perfectly  elastic  manner.  Thus,  the  behavior  of  salt  in  the 
range  of  available  free  field  measurements,  to  strains  as  low  as  approximately 
7  *  10“6,  cannot  be  regarded  as  perfectly  elastic  (Trulio,  1978.  1981). 
Because  of  the  uncertainties  regarding  the  amount  of  energy  that  can  be 
dissipated  as  a  seismic  pulse  propagates  at  intermediate  strains,  Bache,  et 
al ,  (1981)  have  questioned  the  usefulness  of  reduced  displacement  potential 
(RDP)  calculations  based  on  close  range  data  for  the  purpose  of  defining  a 
seismic  source  function. 

The  exact  nature  of  attenuation  in  the  near  field  of  an  explosive 
seismic  source  remains  a  controversial  subject.  In  fact,  the  evidence 
regarding  the  issue  of  linearity  Is  ambiguous.  Using  small-scale  chemical 
explosions  in  pressed  salt  in  the  laboratory,  Larson  (1982)  measured  a  Q  of 
12.5  near  10"3  strain,  and  a  Q  of  24.9  near  6  *  10"4  strain.  The  fact  that  Q 
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is  amplitude  dependent  is  evidence  for  nonlinear  response.  However,  Larson 
also  demonstrated  an  approximately  linear  superposition  of  waveforms  at 
strains  higher  than  10"*,  which  suggests  near-linearity.  The  experimental 
results  reported  in  Tittmann  (1983)  using  resonating  bars  of  natural  halite 
indicate  that  the  linear  anelastic  Q  of  halite  is  quite  high,  near  500,  with  a 
transition  to  an  amplitude  dependent  nonlinear  Q  at  strain  amplitudes  greater 
than  approximately  2  *  10"®.  Burdick,  et  al  (1984a)  have  argued  that  it  is 
possible  to  model  a  seismic  source  function  for  the  Amchitka  tests,  detonated 
in  volcanic  material,  assuming  linear  behavior  in  the  near  field  just  outside 
the  spall  zone  (approx.  700-1200  m/kT^).  Furthermore,  Burdick,  et  al , 
(1984b)  contend  that  the  same  model  can  be  used  to  predict  the  first  vertical 
pulse  arrival  and  rise  times  even  within  the  spall  zone.  They  used  the 
concept  of  a  compressional  elastic  radius  that  in  fact  mey  be  consideraoly 
smaller  than  a  tensional  elastic  radius  and  that  must  extend  at  least  as  far 
as  the  outer  limits  of  the  spall  zone.  Minster  and  Day  (1985)  recently  re¬ 
examined  the  COWBOY  data  set,  and  concluded  that  it  is  possible  to  explain 
simultaneously  the  radial  decay  of  peak  displacement  and  peak  velocity  by 
either  (a)  a  linear  anelastic  model  with  low  Q  (approx.  20)  or  (b)  a  nonlinear 
model  with  amplitude  dependent  Q.  MacCarter  and  Wortman  (1985)  conclude  that 
the  free  field  motion  measurements  from  the  SALMON  event  are  consistent  with 
an  amplitude  independent  Q  of  about  10.  In  any  case,  these  amplitude  inde¬ 
pendent  Q  values  are  too  low  to  represent  the  linear  anelastic  Q  of  salt 
according  to  the  results  of  Tittmann  (1983). 

While  generally  it  is  acknowledged  that  Q  is  defined  rigorously  only 
when  the  material  through  which  a  wave  propagates  behaves  linearly,  it  is 
common  (McCarter  and  Wortman,  1985;  Minster  and  Day,  1985)  to  assume  for 
theoretical  purposes  that  "nonlinear  Q"  can  be  defined  using  the  equations  of 
Mavko  (1979)  and  Stewart,  et  al  (1983): 

Q'1(w,e)  *  Qq”1  ( o>)  +  ae  (1) 
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where  Q"1  is  nonlinear  Q,  Q0-1  is  linear  anelastic  Q,  a  is  a  constant,  and  e 
is  strain.  This  equation  considers  the  combined  effects  of  Coulomb-type 
friction  from  many  intergranular  contacts  on  the  energy  dissipated  during  one 
full  elastic  wave  cycle.  In  their  derivations,  Q”*  was  defined  as 

Q"1  =  (1/2 n) (  £W/W)  (2) 


where  is  the  energy  dissipated  in  one  full  cycle,  and  W  is  the  peak  strain 
energy  stored  per  cycle. 

One  technique  used  in  the  laboratory  to  estimate  nonlinear  0  involves 
the  forced  resonance  of  bars.  As  pointed  out  in  previous  reports  (Tittmann, 
1983a,  1983b,  1984),  the  various  resonating  bar  techniques  offer  a  partic¬ 
ularly  sensitive  tool  for  detecting  the  onset  of  nonlinearity,  even  at  very 
low  strains.  The  data  reported  by  Tittmann  (1983a)  indicate  that  the  intrin¬ 
sic  linear  attenuation  of  dome  salt  is  quite  low  (Q  >  500),  and  that  the  onset 
for  nonlinear  behavior  is  near  2  *  10“6  strain,  corresponding  to  approximately 
1  bar  of  stress.  Similar  behavior  has  been  reported  for  other  rock  types, 
including  granite,  limestones,  sandstones,  and  miscellaneous  igneous  rocks 
(cf  Mavko,  1979;  Stewart  et  al,  1983).  These  measurements  suggest  that 
nonlinear  response  should  persist  to  large  scaled  distances  from  explosions, 
on  the  order  of  104  m/kt*^.  Furthermore,  the  low  Q  values  calculated  by 
Larson  (1982)  [laboratory  measurements],  Trulio  (1979,  1981)  [COWBOY],  Minster 
and  Day  (1985)  [COWBOY]  and  McCarter  and  Wortman  (1985)  [COWBOY,  SALMON]  from 
free  field  measurements  in  salt  are  much  lower.  This  constitutes  additional 
evidence  of  nonlinearity  in  available  free-field  measurements. 

It  is  also  quite  tempting  to  use  the  nonlinear  Q  measurements 
reported  from  forced  resonance-type  experiments  to  estimate  the  coefficient 
in  Eq.  (1).  There  are  certain  dangers  in  doing  this,  however.  First,  it  must 
be  recognized  that  resonating  bar  specimens  are  not  subjected  to  homogeneous 
strain.  This  is  not  a  problem  as  long  as  material  response  is  linear.  How¬ 
ever,  it  tends  to  lead  to  an  underestimation  of  the  value  of  a.  Minster  and 
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Day  (1985)  indicate  a  that  from  flexural  resonance  measurements  is  probably 
underestimated  by  about  a  factor  of  2,  and  from  torsional  measurements  it  is 
underestimated  by  about  a  factor  of  3.  A  second  hazard  with  using  nonlinear  Q 
values  from  resonance  measurements  is  related  to  the  fact  that  Q  is  calculated 
from  the  bandwidth  of  the  resonance  peak,  assuming  linear  anelastic  behav¬ 
ior.  The  distortion  of  the  resonance  peak  shape  associated  with  nonlinear 
resonances  has  been  described  by  Tittmann  (1983b)  and  casts  an  element  of 
suspicion  on  the  significance  of  the  nonlinear  Q  values  reported.  Finally, 
during  both  extensional  and  flexural  resonances  each  increment  of  volume 
within  the  specimen  is  subjected  to  tensile  stresses  half  of  the  time  and  to 
compressive  stresses  the  other  half.  A  tacit  assumption  is  that  the  material 
behaves  the  same  in  compression  as  it  does  in  tension.  The  apparent  drop  in 
bar  modulus  with  increasing  vibration  amplitude  has  been  noted  (Tittmann, 
1983a,  1983b,  1984),  and  is  probably  also  related  to  the  distortion  of  the 
resonance  peak  shape  with  increasing  strain.  Not  proven,  however,  is  whether 
either  the  frictional  losses  or  the  effective  modulus  of  a  rock  are  the  same 
in  tension  as  they  are  in  compression.  Resonating  bar-type  measurements  are 
not  capable  of  resolving  the  details  of  the  rheological  response  of  materials 
to  nonlinear  loading. 

On  the  basis  of  the  above  discussion,  the  laboratory  work  most 
beneficial  to  resolving  the  issue  of  linearity  vs  nonlinearity  would  involve 
studies  of  the  mechanical  response  of  various  test  site  materials,  including 
salt,  to  realistic  seismic  pulse  loading  histories.  Such  studies  would  neces¬ 
sarily  require  the  application  of  true  uniaxial  strain  while  test  specimens 
are  subjected  simutaneously  to  a  hydrostatic  bias  stress.  This  type  of  study 
is  very  difficult  technically,  and  results  of  such  studies  do  not  exist,  at 
least  not  in  the  open  literature.  Several  papers  have  examined  the  details  of 
stress-strain  hysteresis  loops  under  nonlinear,  uniaxial  stress  loading  condi¬ 
tions  (distinguished  from  uniaxial  strain)  (Gordon  and  Davis,  1968;  McKavanagh 
and  Stacey,  1974),  but  in  each  of  these  cases  only  congressional  measurements 
were  performed.  Shock  and  Louis  (1982)  examined  the  inelastic  response  of  two 
rock  types,  sandstone  and  granite,  to  high  tensile  stress  loading,  but  did  not 
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examine  hysteresis,  per  se.  Nor  was  the  behavior  in  axial  compression  con¬ 
trasted  with  the  behavior  in  tension. 

In  this  report  we  describe  the  results  of  combined  high  amplitude 
tensile  and  compressive  uniaxial  stress  loading  experiments  on  four  different 
rock  types:  Westerly  granite,  Boise  sandstone,  Berea  sandstone,  and  Indiana 
limestone,  under  ambient  pressure  conditions.  This  represents  the  first  stage 
of  a  new  effort  with  the  ultimate  goal  to  examine  the  mechanical  response  of 
test  site  materials,  including  dome  salt,  to  high  amplitude  uniaxial  strain 
loading  at  elevated  hydrostatic  pressures.  We  examine  the  details  of  the 
stress-strain  hysteresis  loops,  and  demonstrate  that  both  the  elastic  and  the 
inelastic  properties  of  rocks  at  nonlinear  amplitudes  are  quite  different  in 
compression  than  in  tension.  Results  to  date  indicate  that  the  behavior  in 
compression  at  strains  greater  than  10~5  is  nonlinear,  but  with  much  less 
hysteresis  than  in  tension.  The  behavior  in  tension  at  comparable  strains  is 
also  nonlinear,  but  a  relatively  large  hysteresis  loop  is  evident.  This 
indicates  that  frictional  losses  are  much  greater  in  tensile  loading  than  in 
compressive  loading. 
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3.0  EXPERIMENTAL  PROCEDURES 

Figure  3.1  includes  schematic  illustrations  of  the  mechanical  parts 
of  the  system  used  for  measurement.  Briefly,  cylindrical  rock  specimens 
2.54  cm  in  dia.  and  6.35  cm  long  Mere  cut  from  blocks  of  Boise  sandstone, 
Berea  sandstone.  Westerly  granite,  and  Indiana  limestone.  A  resistance-type 
strain  gauge  was  bonded  to  the  sur.face  of  each  specimen  to  respond  to  strains 
in  the  longitudinal  direction.  The  rock  was  then  bonded  adhesively  to  stain¬ 
less  steel  mounting  fixtures  using  a  semiplastic  polymer  resin  (Crystal  Bond 
#509,  Aremco  Products)  that  is  quite  brittle  and  stiff  at  room  temperature, 
but  that  softens  readily  at  approximately  70°C.  By  keeping  the  specimen  warm 
while  mounting  it  in  the  apparatus,  residual  stresses  could  be  removed. 
Measurements  were  made  at  room  temperature  and  under  ambient  pressure 
conditions. 


MTS  LOAD  FRAME 


SAMPLE 


STRAIN  GAUGE 


ROCK  SPECIMEN 


EXTENSION  BARS 


ADHESIVE  BONDING 


Fig.  3.1  Schematic  illustration  of  the  apparatus  used  for  the  measurement 
of  stress-strain  response  curves. 
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Experiments  were  run  using  an  MTS  electro-hydraulic,  closed-loop  load 
frame  equipped  with  the  hydraulic  hardware  and  controls  necessary  to  apply 
alternating  tensile  and  compressive  loads.  The  instrumentation  used  for  the 
measurements  is  illustrated  in  Fig.  3.2.  The  frequency  of  loading  was  usually 
1  Hz,  except  for  two  series  of  measurements  on  Westerly  granite  and  Berea 
sandstone  at  0.1  Hz.  A  series  of  eight  cycle  bursts  of  load-controlled  con¬ 
stant  strain  rate  triangle  wave  loading  were  applied  to  each  specimen. 
Measurements  of  load  and  strain  were  collected  simultaneously,  digitized,  and 
stored  in  the  computer. 

Representative  time  functions  of  stress  and  strain  are  illustrated  in 
Fig.  3.3.  Measurements  were  obtained  on  each  sample  using  bursts  of  increas¬ 
ing  maximum  stress,  starting  with  t  3.4  *  105  Pa  and  working  up  to  a  maximum 
of  +  3.4  *  10^  Pa,  unless  precluded  by  the  breakage  of  the  specimen.  As  many 
as  5  bursts  of  8  waves  each  were  applied  at  low  amplitudes  in  order  to  enable 
satisfactory  signal  averaging  when  the  signal -to-noise  ratio  was  low.  At 
higher  amplitudes  the  signal -to-noise  ratio  was  good,  and  only  one  eight-cycle 
burst  was  run  in  order  to  avoid  excessive  damage  to  the  specimen. 
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Fig.  3.2  Schematic  illustration  of  Instrumentation  used  to  acquire  and 
process  stress-strain  response  curves. 
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Fly.  3.3  Representative  stress  loading  functions  with  corresponding  strain 
response  curves,  (a)  MTS  electro-hydraulic  load  fram  and  (b)  detail 
of  sample  and  end  fixtures. 


10 

C7069TC/bje 


Rockwell  International 

Sc  tone*  Center 


SC5361.8SA 


4.0  EXPERIMENTAL  RESULTS 

To  test  the  linearity  of  the  experimental  apparatus  and  the  stiffness 
of  the  adhesive  bonds,  the  first  phase  of  the  experimental  study  involved 
measurements  of  stress  and  strain  on  an  aluminum  bar  which  has  a  relatively 
ideal  linear  elastic  response.  A  typical  curve  of  stress  vs.  strain  for  an 
aluminum  bar  is  shown  in  Fig.  4.1.  Extensional  stresses  and  strains  are 
positive.  Only  a  small  amount  of  hysteresis  is  observed  in  this  set  of 
measurements,  and  the  apparent  relationship  between  stress  and  strain  is 
nearly  linear.  A  small  amount  of  hysteresis  is  observed  and  may  be  attributed 
to  plastic  flow  in  the  resin  used  to  bond  the  aluminum  bar  to  the  loading 
fixtures.  This  is  a  relatively  insignificant  effect  compared  with  the  large 
amounts  of  hysteresis  and  nonlinearity  observed  in  rocks  under  corresponding 
stresses. 


SRMPLF  Sljormcc  STPPIN 

Fig.  4.1  Stress-strain  response  for  aluminum  test  bar. 

Experimental  results  also  show  that  each  rock  displays  nearly  linear 
behavior  when  the  loading  is  very  small,  around  4.5  *  104  Pa,  resulting  in 
strains  between  2  *  10”®  and  6  *  10”®,  depending  on  rock  type  (Figs.  4.2 
through  4.5).  At  higher  load  levels  (Figs.  4. 6-4. 9)  there  is  evidence  of 
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strong  nonlinearity  and  inelasticity  in  the  weakest  frame  rock,  Berea  sand¬ 
stone,  while  the  nonlinear  and  inelastic  effects  are  more  subtle  in  the 
stronger  frame  rocks,  especially  Indiana  limestone.  In  all  rock  specimens 
except  Indiana  limestone  the  effective  modulus  of  the  rock  is  significantly 
larger  in  compression  than  in  tension.  In  two  anomalous  cases,  Indiana 
limestone  and  Westerly  granite  have  a  slightly  higher  modulus  in  tension  than 
in  compression  at  rather  small  loads  (Fig.  4.10). 

McKavanagh  and  Stacey  (1974)  argue  that  it  is  possible  to  distinguish 
between  linear  anelastic  relaxation  and  nonlinear  relaxation,  such  as  inter¬ 
granular  friction,  by  examining  the  shape  of  the  hysteresis  loop  tips.  Cusped 
tips  indicate  a  nonlinear  mechanism,  while  rounded  tips  indicate  linear 
anelastic  relaxation.  Expanded  views  of  the  hysteresis  loop  tips  measured  on 
Westerly  granite,  Boise  sandstone,  and  Berea  sandstone  are  shown  in  Figs. 

4.11,  4.12,  and  4.13,  respectively.  In  each  case  the  tips  are  cusped  and  not 
rounded,  which  argues  against  an  anelastic  relaxation  mechanism. 

In  some  cases  evidence  for  large  nonrecoverable  changes  in  the  rock 
during  the  first  excursion  to  a  "new"  maximum  load  is  clear.  A  reasonable 
assumption  is  that  these  changes  are  associated  with  intergranular  cracking. 
This  is  apparent  when  the  hysteresis  loop  is  not  clearly  established  until  at 
least  the  second  cycle  of  loading.  In  the  case  of  Berea  sandstone  subjected 
to  1.25  *  105  Pa  axial  stress  (Fig.  4.14),  the  result  of  the  first  full 
loading  cycle  was  a  net  shortening  of  the  specimen.  The  nonrecoverable  change 
at  comparable  stress  is  small  in  rocks  with  a  moderately  strong  frame,  such  as 
Boise  sandstone  and  Westerly  granite  (Figs  4.15  and  4.16).  Indiana  limestone, 
the  rock  with  the  strongest  frame,  shows  no  hysteresis  and  no  nonrecoverable 
changes  after  comparable  loading  (Fig  4.17). 
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BOISE  SANDSTONE  B60-34B 


Fig.  4.2  Stress-strain  response  for  Boise  sandstone  with  a  peak  loading 
stress  of  4.5  *  10 **  Pa  showing  nearly  linear  behavior. 


BEREA  SANDSTONE  RS02-20B 


SAMPLE  SURFACE  STRAIN 

Fig,  4.3  Stress-strain  response  for  Berea  sandstone  with  a  peak  loading 
stress  of  4.5  *  10**  Pa  showing  nearly  linear  behavior. 
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Fig.  4.4  Stress-strain  response  for  Westerly  granite  with  a  peak  loading 
stress  4.5  *  101*  Pa  showing  nearly  linear  behavior. 


INDIANA  LIMESTONE  333-4 


SAMPLE  SURFACE  STRAIN 


Fig.  4.5  Stress-strain  response  for  Indiana  limestone  with  a  peak  loading 
stress  of  4.5  *  10 **  Pa  showing  nearly  linear  behavior. 
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Fig.  4.6  Stress-strain  response  for  Boise  sandstone  with  a  peak  loading 
stress  of  4.5  *  10H  Pa. 
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Fig.  4.7  Stress-strain  response  for  Berea  sandstone  with  a  peak  loading 
stress  of  1.32  *  10 6  Pa. 
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Fig.  4.8  Stress-strain  response  for  Westerly  granite  with  a  peak  loading 
stress  of  1.32  *  10 6  Pa. 
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Fig.  4.9  Stress-strain  response  for  Indiana  limestone  with  a  peak  loading 
stress  of  1.32  *  10 6  Pa. 
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Fig.  4.10  Stress-strain  response  for  two  rocks  showing  the  anomalous 
occurrence  of  a  slightly  higher  effective  modulus  in  tension 
than  In  compression.  Peak  loading  stresses  are  relatively 
small  (1.32  *  10s  Pa). 
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.11  Magnified  views  of  the  hysteresis  loop  tips  for  Westerly  granite. 
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4.12  Magnified  views  of  the  hysteresis  loop  tips  for  Boise  sandstone. 
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Fig.  4.13  Magnified  views  of  the  hysteresis  loop  tips  for  Berea  sandstone. 
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SRMPLE  SURFRCE  STRRIN 

Fig.  4.14  Comparison  of  the  first  loading  cycle  for  Berea  sandstone  at 

1.25  *  10  6  Pa  with  the  subsequent  7  cycles  showing  nonrecoverable 
changes  after  the  first.  The  previous  maximum  load  was 
4.39  x  10 5  Pa. 
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BOISE  SRNDSTONE  B80-34R 
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Fig.  4.15  Stress-strain  response  for  the  first  eight  loading  cycles  for 

Boise  sandstone  at  1.25  *  10 b  Pa.  A  small  amount  of  nonrecover- 
able  change  is  observed  after  the  first  cycle,  but  the  following 
seven  are  superimposed.  The  previous  maximum  load  was 
4.39  x  10s  Pa. 
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Fig.  4.16  Stress-strain  response  for  the  first  eight  loading  cycles  for 

Westerly  granite  at  1.25  *  10 6  Pa.  A  small  amount  of  nonrecover- 
able  change  Is  observed  after  the  first  cycle,  but  the  following 
seven  are  superimposed.  The  previous  maxlmun  load  was 
4.39  *  10a  Pa. 
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Fig.  4.17  Stress-strain  response  for  the  first  eight  loading  cycles  for 
Indiana  limestone  at  1.2b  *  10b  Pa.  The  material  appears  to  be 
nearly  linear  elastic  and  nonrecoverable  changes  are  not  observed. 

By  reversing  the  direction  of  initial  loading,  the  shape  of  tensional 
hysteresis  loops  and  compressive  hysteresis  loops  can  be  examined.  In  Figs. 
4.18  and  4.21  we  show  the  effects  of  reversing  the  loading  direction  from 
tension-first  to  compression-first,  and  from  compression-first  to  tension- 
first,  respectively,  for  two  different  rock  types.  It  is  apparent  that  in 
Boise  sandstone  the  opening  of  the  loop  is  due  primarily  to  tensile  stresses, 
and  that  the  behavior  in  compression  is  nearly  nonlinear  elastic.  In  Westerly 
granite  these  effects  are  similar,  but  the  contrast  between  compressional 
hysteresis  and  tensional  hysteresis  is  less  pronounced.  A  hysteresis  loop 
exists  in  both  compression  and  tension,  but  it  is  smaller  in  compression. 

Finally,  the  effect  of  frequency  on  the  shape  of  the  hysteresis  loop 
has  been  examined.  The  results  of  measurements  on  Boise  sandstone  and  Wes¬ 
terly  granite  at  1  Hz  and  at  0.1  Hz  are  shown  In  Figs.  4.22  and  4.23.  The 
curves  superimpose  very  well,  which  Is  consistent  with  a  frictional  relaxation 
mechanism,  even  though  this  observation  alone  does  not  preclude  anelastic 
relaxation. 
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Fig.  4.18  Effect  of  reversing  the  Initial  loading  direction  from  tension- 
first  to  compression-first  for  Boise  sandstone  with  a  maximum 
loading  stress  of  4.39  *  10b  Pa,  showing  the  lack  of  significant 
hysteresis  In  the  compresslonal  part  of  the  cycle. 
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Fig.  4.19  Effect  of  reversing  the  Initial  loading  direction  from  compression 
first  to  tension  first  for  Boise  sandstone  with  a  maximum  loading 
stress  of  4.39  *  10 6  Pa,  showing  the  development  of  a  hysteresis 
loop  In  tension. 
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Fig.  4.20  Effect  of  reversing  the  Initial  loading  direction  from  tension- 
first  to  compression-first  for  Westerly  granite  with  a  maximum 
loading  stress  of  4.39  *  10®  Pa,  showing  a  significant  congres¬ 
sional  hysteresis  loop. 
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Fig.  4.21  Effect  of  reversing  the  Initial  loading  direction  from  compression 
first  to  tension-first  for  Westerly  granite  with  a  maximum  loading 
stress  of  4.39  *  10 6  Pa.  Figure  shows  the  existence  of  a  hystere¬ 
sis  loop  In  tension  which  Is  larger  than  the  compresslonal  hystere 
sis  loop  Illustrated  In  Fig.  4.20. 
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Fig.  4.22  Comparison  of  the  stress-strain  response  for  Berea  sandstone  at 
two  different  frequencies.  The  traces  appear  to  be  Identical. 
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Fig.  4.23  Comparison  of  the  stress-strain  response  for  Westerly  granite  at 
two  difference  frequencies.  The  traces  appear  to  be  Identical. 
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5.0  DISCUSSION  AND  CONCLUSIONS 

In  this  study  we  have  examined  the  details  of  stress-strain  hystere¬ 
sis  loops  when  rock  specimens  are  subjected  to  alternating  compressive  and 
tensile  stresses.  In  general  at  nonlinear  amplitudes  the  stiffness  of  the 
rock  is  greater  in  compression  than  in  tension.  Furthermore,  preliminary 
results  indicate  that  most  of  the  energy  loss  during  a  full  cycle  of  loading 
occurs  as  a  result  of  strain  in  extension;  the  hysteresis  loop  in  compression 
is  smaller  than  the  hysteresis  loop  in  tension.  The  shape  of  the  hysteresis 
loops  also  appears  to  be  independent  of  frequency.  These  observations  indi¬ 
cate  a  loss  mechanism  associated  with  intergranular  friction.  Intergranular 
sliding  appears  to  be  restricted  by  the  impingement  of  opposing  crack  faces  in 
compression,  since  the  rock  is  stiffer  in  compression  than  in  tension,  and  a 
large  hysteresis  loop  develops  only  when  the  rock  is  subjected  to  tensile 
stresses. 

The  results  of  this  study  indicate  that  the  mechanical  behavior  of 
rocks  can  be  significantly  different  in  compression  than  in  tension  and  that 
the  onset  of  nonlinear  effects  with  increasing  strain  may  not  be  the  same  for 
tensile  loads  as  for  compressive  loads.  All  available  evidence  indicates  that 
the  primary  relaxation  mechanism  at  nonlinear  amplitudes  between  10"®  strain 
and  10"*  strain  Involves  intergranular  friction.  More  experimental  work  in 
this  area  will  shed  light  on  the  Issue  of  linearity  vs  nonlinearity  at  inter¬ 
mediate  strains,  and  also  will  provide  realistic  detailed  Information  about 
rock  rheology  for  the  numerical  modeling  of  near-field  seismic  pulse  propa¬ 
gation.  For  this  work  to  be  most  meaningful  true  uniaxial  strain  is  neces¬ 
sary,  and  specimens  must  be  exposed  to  elevated  confining  pressures. 
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